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We present a method for constructing stationary, asymptotically flat, rotating 
solutions of Einstein's field equations. One of the spun-up solutions has quasilocal 
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I mass but no global mass. It has an ergosphere but no event horizon. The angular 



q ' momentum is constant everywhere beyond the ergosphere. The energy-momentum 

content of this solution can be interpreted as a rotating string-fluid. 
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q-i A metric description of the spacetime around an object encodes information about its 



I. INTRODUCTION 



mass and angular momentum. While many of the objects of interest to astrophysicists 
and cosmologists are rotating, there are not a large number of metrics that can be used 
to describe isolated rotating objects. Most models of the early universe encode no angular 

b : 

momentum into the metric. Developing methods to generate rotating solutions is clearly of 
value. There are two well known solution generating algorithms which start with a static 
metric and transform it by adding rotation to some aspect of the spacetime structure. The 
complex transformation discovered by Newman (N-J) and others [lj], for example, 

takes the Schwarzschild metric to the Kerr metric, and the Reissner-Nordstrom metric to 
the Kerr-Newman solution. The N-J transform adds angular momentum to a seed spacetime 
initially containing only global mass (and possibly charge). The N-J method has also been 
used |4j to spin up spherical metrics to obtain Kerr interiors. There is a method begun 
by Ehlers Jij] and fully developed by Geroch which adds twist to the timelike Killing 
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vector of the original (seed) metric and also changes its norm. If the method is applied to 
the static Schwarzschild spacetime the Taub-Nut spacetime results. 

In this paper, we suggest a new method for adding rotation to a known static spacetime. 
The method starts with an asymptotically flat static metric written in null Bondi-Sachs 
form. A transform of the type dip = dip — Vtdu is applied. Q has coordinate dependence thus 
adding angular momentum to the spacetime. The Komar integral is used to calculate the 
angular momentum. Asymptotic flatness is enforced by using the Bondi-Sachs metric as a 
prototype for adding angular momentum, with the Bondi-Sachs metric conditions ensuring 
that asymptotic flatness is preserved in the transformation. 

In the next section we describe the method in detail. In the third part of the paper we 
apply the method to the Schwarzschild spacetime, adding a rotation term, maintaining the 
Bondi-Sachs metric form, and using the rotation term from the linearized Kerr metric in the 
Bondi-Sachs frame as a guide. We also apply the method to Minkowski space, obtainable 
as a limiting case of the spinning Schwarzschild example. This spacetime contains some 
unusual features. Unusual because there is quasilocal mass but no global mass, and there 
is an ergosphere but no event horizon. The solution is stationary, axially symmetric, and 
asymptotically flat. The energy-momentum content can be identified as a rotating string 
fluid. Values from the Komar superpotential provide angular momentum which does not 
fall off asymptotically and quasilocal mass which does. The quasilocal mass arises from the 
kinetic energy of the rotating string fluid. 

CONVENTIONS 

In this work Greek indices range over (0,1,2,3) = (u,r,$,(p). Sign conventions are 
2A u .[ a/ 3] = ApiR^ap, and R a p = R u a/3u - The metric signature is (+,-,-,-) and the field equa- 
tions are G, w = —8ttT uv . 



II. THE METHOD 



The Bondi-Sachs metric has been used since the 1960s. It is based on a foliation of 
outgoing null hypersurfaces (A/") and a conformal boundary (X), with sufficient generality 
to describe bounded radiating astrophysical systems (details are given in Appendix D). 
Asymptotically flat electrovac spacetimes in Bondi-Sachs form are summarized by Bicak 
and Pravdova |8f. 
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To construct our solution, we begin with the Schwarzschild vacuum metric. Angular 
momentum is added using the Killing vector as generator of rotations. This is done 
in the context of the Bondi-Sachs metric in order to identify the term to add, and to 
maintain asymptotic flatness. A new term is added containing angular velocity with form 
dip = d(p — Vtdu. For Q a function rather than a constant, this is not an integrable coordinate 
transformation. The Komar superpotential integral for d v insures that the desired amount 
of angular momentum is added. The particular choice of fl is guided by the linearized Kerr 
metric written in Bondi-Sachs form. We then have an analytic asymptotically flat spinning 
metric. 



A. Kerr Vacuum Solution 

The Kerr solution in Boyer-Lindquist coordinates is given by 

glydx a dx p = ip dt 2 - (E/A) dr 2 + (1 - V)2asin 2 tf dtdip 
— £ d$ 2 - sin 2 tf[£ + (2 - ^)a 2 sin¥] d V 2 

Here if) = I — 2m r/S, E = r 2 + a 2 cos 2, $, A = r 2 + a 2 — 2m r. The linearized Kerr metric 
has terms linear in 'a'. 

glj kerr dx a dxP = (1 - 2m /r)dt 2 - (1 - 2m Q / r)- 1 dr 2 + (2m Q a/r)sm 2 $ 2dtd V (1) 
-r 2 (^ 2 + sin 2 ^)V- 

To reach Bondi-Sachs form transform from Boyer-Lindquist to outgoing u = t — r — 
2m ln(r — 2m ). 

g l ^- kcrr - BS dx a dx 13 = (1 - 2m /r)du 2 + 2dudr + (2m a/r)sin 2 tf 2dudtp (2) 
-(1 - 2m /ry 1 (2m a/r)sm 2 $ 2drdip - r 2 (dtf 2 + sin 2 'd) dip 2 . 

Note that 

{du) a {du)p g^ TI = to order V. 
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III. SPINNING SCHWARZSCHILD 

The Schwarzschild vacuum is spun up by choosing a term with angular velocity Q(r). 
The spinning metric is written in Bondi-Sachs form as 

g^dx a dx 13 = (1 - 2m /r)du 2 + 2dudr - rW - r 2 sin 2 (dip - Qdu) 2 . (3) 

The term (2m a/r)sin 2 $ 2dudip in the linearized Kerr metric Eq.(J2J) leads to the function 
U v = n(r) in Eq.(jDU), multiplied by r 2 E w 

— r 2 g j n ^ This directs the choice Q(r) = r^/r 3 . 
The Bondi-Sachs form can be rewritten as the spinning Schwarzschild (sSCH) metric 

gf^ K dx a dx p = Zdu 2 + 2dudr + 2fir 2 sin 2 tf dudip - r 2 (dd 2 + sin 2 d dip 2 ) (4) 

with Z :— 1 — 2mo/r — (f2rsin#) 2 . 

To understand the physical content of the sSCH solution we write the metric in a basis 
which is locally non- rotating (v a has zero twist). 

gfp H = v a V/3 ~ r a rp - ftcftp - Cpjpp. (5) 

Here A 2 = 1 — 2m /r. The unit vectors are defined by 

v a dx a = A~ 1 (A 2 du + dr), v a d a = A^ 1 (d u + nd v ), (6a) 

r a dx a = A- X dr, f a d a = A~ x (d u - A 2 d r + £ld v ), (6b) 

$ a dx a = rdd, d d a = —r~ 1 d$, (6c) 

& a dx a = rsm$(dip — Vtdu), (p a d a = — (rshrt?) -1 ^. (6d) 

The Einstein tensor is expanded as 

Gf™ = U 2 {v a v p - r a r p + + Z(p a {p p ) (7) 

with U 2 = (9/4)n 2 sin 2 tf. In the locally non-rotating frame, the energy-momentum tensor 
is 

T ap CR = pVaVp + p r r a rf3 + Pv&a'dp + Vtpfyjpp- ( 8 ) 
It follows from G a p = —8>-nT a p that 

8np = -U 2 , 87r Pr = U 2 , (9) 
8irp# = -U 2 , 8-k Pv = -3U 2 . 
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The Ricci tensor has the form 

Rf™ = 2U 2 (v a Vp-r a rp + <p a <pp), R% cn f = 2U%, (10) 

with l a d a = d r an eigenvector of the Ricci tensor. 

The fluid velocity v a is expansion- free and twist-free, with shear scalar \a a pa al3 = U 2 , 
and acceleration a a = (-^)r a . The content of spinning Schwarzschild can be interpreted as 
a rotating radial string fluid with p + p r = 0. 

The sSCH metric can be written as 

sS^-sSF-WaSp, (11) 

where is the Schwarzschild metric and Sp is a spacelike vector such that 

Spdx 13 = (VL/2)r 2 sin 2 Mu-r 2 sin 2 &dtp, (12) 
S a S a = -(rsintf) 2 . 

When m — > in sSCH, we have spinning Minkowski (sM): 

C = V«p - 20l a Sp. (13) 
This spacetime, as seen from null infinity, has angular momentum but no mass. 

A. Spinning Minkowski 

We set m = in Z of Eq.(J3J) to obtain a metric with no global mass. The spinning 
Minkowski metric (sM) has one-parameter, t~o, in function Q(r) = r 2 /r 3 . 

gfp l dx a dx p = (1 - Q 2 r 2 sin 2 $)du 2 + 2dudr + 2fir 2 sin 2 tf dud<p - r 2 {d$ 2 + sin 2 $ dtp 2 ). (14) 

The metric has two Killing vectors, stationary kr n ) and axial kr^y In terms of the Minkowski 
metric rj a p, the sM metric can be written as 

g s ^dx a dx p = T] a pdx a dx> 3 - tt 2 r 2 sin 2 $du 2 + 2ttr 2 sin 2 $ dud<p. 

The rotation axis, 0, 2tt] is a regular line in the spacetime. As r — > oo, Qr 2 — > 
and g^p — > r] a/3 , hence the sM metric is asymptotically flat. In fact, one can remap the 
Minkowski metric 

rj a pdx a dx p = du 2 + 2dudr - r 2 (d$ 2 + sin 2 $d& 2 ) 
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by shifting the p coordinate to p = p — ttoii. 

r ]ap dx a dx 13 = (1 - n 2 r 2 sm 2 $)du 2 + 2dudr + 2fi r 2 sin 2 7? dudip - r 2 (dtf 2 + sin 2 dp 2 ). 

The Minkowski metric then coincides with the limit of the sM metric as ^(J, 1 ™^) — > 0, and 
tt -> 0. 

The sM solution is written in the same locally non-rotating frame as sSCH above, but 
now the frame velocity v a is geodesic. It is also expansion-free and twist-free, with shear 
scalar \o a ^ = U 2 = (3r 2 sin^/2r 3 ) 2 . 

& = Vcfip - r a r p - tijtp - p a p p . (15) 

The unit vectors are [this is the tetrad in Eq.© with A = 1] 

v a dx a = du + dr, v a d a = d u + Qd v , (16a) 

f a dx a = dr, f a d a = d u — d r + Qd^, (16b) 

# a dx a = rdtf, $ a d a = -r~ l dv, (16c) 

p a dx a = rsnxd{dp> — Vtdu), p a d a = —(rsmd)~ 1 d ip . (16d) 

The Einstein tensor and energy- momentum tensor here are the same as those of the sSCH 
solution above. 

G% = U 2 {v a vp - f a f p + 0j p + ZCp a Cpp). (17) 

T af = PVaV/3 + PrfaTfi + P^a^/3 + p v p a p^ (18) 

8np = -U 2 , Sixp r = U 2 , (19) 
8tt^ = -U 2 , &irp v = -3U 2 . 

The sM solution contains a rotating radial string fluid with p + p r = 0. Since U = 
(3/2)ro sin^/r 3 , the entire matter content approaches vacuum as r — > 00. 

IV. MASS 

The Bondi mass of g sSCH is evaluated on a spherical cut of future null infinity X + . Equa- 
tion ()D4j) and the Weyl tensor component Eq. (|A2J) yields 

M bondl = ~ <f(*°2 + *°) >/=gdMp = m . 

dAf 
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Since metric g sM is the m — > limit of g sSCH it follows that 

MSL = 0- (20) 

The Komar superpotential for Killing vector Ar is 

U aP {k) = (- 9 )5[V¥ - V p k a \. (21) 

For timelike Killing vector k^sdp = d u one writes U a,3 (d u ), and for axial Killing vector 
k^dp = d v one writes U a ^(d (p ). Metrics g sSCH and g sM both have (— g)z = r 2 sm$. For 
u = const three-surface M, the mass within a u = const, r = const two-surface DM is 

M komar = ~<f U a ' 3 (d u ) dS a p. (22) 



in 
av 



Metric g^ CH and two-surface dS a p = u,[ a r,p} dfidip contains global mass m and quasilocal 
mass at all r beyond the ergosphere. 



M sSCH 
komar 



27V 7T 

-— / /"(-2m sin^-3^sin 3 ^)c/^ = m + ^4, (23) 
87r J J r 6 r 6 



o o 

Mt M = *X 

komar 3 



There is no global mass for g" M since m = and M^ mar — > as r — > oo. 

Although the sM solution is not spherically symmetric, we can compute a sectional cur- 
vature mass since tetrad vectors d and (p a of Eq. (J16c)) and (|16d|) form a bivector which 
satisfies the Frobenius suface-forming condition. They span a family of closed two-surfaces. 
The sectional curvature (Gaussian curvature) mass is 

2M sSCH 



curv 



R a p^ <p p $y v = mo. (24) 



The sectional curvature mass for g sM vanishes. 

The explanation of why M sSCH and M sM are positive when the density p = —U 2 /87T is 
negative lies with the relativistic contribution of pressure to mass. If one traces the Komar 
integral back it then becomes clear. 

M komar = ju a P{d u )ds a p 



8N 

1 

8^ 

av 

1 

8tt 

Af 
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for dS/3 = u^d 3 x. The three- volume integration is over any u = const null surface from 
(beyond the ergosurface) to oo. The Ricci tensor provides 



Lkomar 



Mi,,,,,,,,,- — /' sf-g{ '2f> l>< I p v )d 3 x 
ovr J 

gU 2 drd'ddip 



1 



Air 

For g sScn the volume integral does not measure the m contribution. 



V. ANGULAR MOMENTUM 

For u = const three-surface Af, the angular momentum within a u = const, r = const 
two-surface dAf is 

J=-^ju a P(d v )dS a p. (25) 

Metric g sSCB - and metric g sM admit constant J at all r beyond the ergosphere. For two-surface 
element dS a/ 3 = u,[ a r,/3] dfidip 

tsSCH tsM 1 / / /r. 2 • 3 „a\ j„q J,„ _ 1 2 



jsbCH = JS M = _ / / (3r 2 gin 3 ^ Mdip = __ r 2 (26) 

167T ./ ./ 2 







This is a global value. 

The notation of the Ko m a, « is chosen Q so that 

J kcrr («9 v ) = -m a and M kerr (<9 n ) = m . 

The sign difference occurs because our —2 signature gives timelike vectors positive norms 
and spacelike vectors negative norms. 

VI. SURFACES OF SM 
A. Horizon 



Spacelike closed two-surfaces are trapped if their null generators (incoming and outgoing 
null geodesies) converge toward the future. The time evolution of the outermost trapped 
surface is a null surface called the " apparent horizon" . The spherical two-surfaces of metric 
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g sM have area 4irr 2 , with null geodesic generators l a and n a [see Eq. ()A4|) ]. Their expansions 
(convergences) are given respectively by 

p = -l a d a (\nr) = -1/r, (27) 
H = n a d a (\nr) = -l/(2r). 

(here p is a spin-coefficient, not mass-density). There can be no horizon since neither p nor 
p change sign over the interval < r < oo. 



B. Ergosurface 

The norm of the timelike Killing vector is, with Q = r^/r 3 , 

d^-df 1 = l-ft 2 r 2 sin 2 tf, (28) 

placing the ergosphere at r 2 = ±r 2 smity =\ sin^ |. 

The shape of the ergosurface (outer boundary of the ergosphere) is determined by the 
locus of points where the norm of <9^ M is zero, r 2 =| Tq sin$ |. The curve ranges from r = 
to r = r , and is called the "Leminscate of Bernoulli" [lpj . In (x,y,z) coordinates, with 
rotation about the z-axis, the ergosurface is shaped like a deformed torus without a center 
hole, since the surface boundary curve goes smoothly to at i = ?/ = 2 = 0. Inside the 
ergosurface the spacelike two-surfaces are two-spheres. 

With <9^ M • <9^ M < 0, the Gaussian curvature of two-surfaces inside the ergosphere is 

/C ergo = 1/r 2 . (29) 

We compute the Euler-Poincare characteristic x °f the ergosurface. The two-surface metric 
is 

dsj_ smf = r 2 sin d(dd 2 + sin 2 d d V 2 ) (30) 

and 

^2-ergo = n • q (31) 

Tq sm v 



undefined at the poles $ = 0, 7r. The Euler-Poincare characteristic 



lll | is then 



2-7T 7T 



2n X = J J /C 2 -cr g o(ro sin 2 #) = 4vr. (32) 
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Thus x — 2, with associated genus g = which characterizes the topology of a standard 
two- sphere. 

The Kerr ergosphere is quite different. The norm of the Kerr timelike Killing vector (in 
Boyer-Lindquist coordinates) is 

^Kerr _ ^Kerr _ ^ 2m r 

1 1 r 2 + a 2 cos 2 1? ' 

The Killing vector is spacelike in the region 2m < r < m + (m 2 , — a 2 cos 2 $) 1//2 . The Kerr 
ergosurface is tangent to the rotating trapped surface at the poles. 



C. Geodesic Deviation of sM 



Test particles move along a congruence of timelike geodesies v a = dx a (r)/dr. Vector v a is 
twist-free and acceleration-free with respect to metric g sM and is given by v a d a = d u + Qd^. 
Covariant differentiation along the geodesies is defined by 

D/dr := v a V a . (33) 

Deviation vector i] a is tangent to the line connecting pairs of neighboring geodesies in the 
congruence. It is Lie transported along the congruence and satisfies 

C vV a = 0. (34) 

jp2i] a gives the relative acceleration between neighboring test particles. The geodesic devi- 
ation equation is 

^ a =[i^tA>V. (35) 

If one considers any ip = const slice of a family of ip} two-spheres then 77 = (1/ sintyd^ is 
a deviation vector linking timelike geodesies orthogonal to the two-spheres. We find 

^ = (hfrW (36) 



for Q = Tq/t 3 and tetrad vector (j!6d|) . A view of the geodesic congruence shows the tip of 
the connecting vector rj spiraling up the congruence with acceleration Q 2 r. 
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VII. DISCUSSION 

A method has been presented for adding angular momentum to static spacetimes. Metric 
gsSCH j g a S pi nri i n g generalization of the vacuum Schwarzschild metric. The spinup method 
creates a spinning object with an atmosphere in an asymptotically flat spacetime. There 
are many isolated astrophysical systems in our observable universe that spin. We have 
recently seen that the unexpected behavior of the angular velocity of globular clusters was a 
strong clue to the existence of dark matter and dark energy. Analytic metrics, with a family 
of possible angular velocities, can be a useful tool for studying and classifying rotating 
systems. They can be used to predict and analyze observable data such as lensing and to 
check numerical simulations. 

Here we started with the Schwarzschild metric and added a rotation term, maintaining 
Bondi-Sachs metric form and using the rotation term from the linearized Kerr metric (trans- 
formed to a Bondi-Sachs frame) as a guide. The resulting metric, g sSCH ; has global mass 
and angular momentum. The global mass is M^pJ 1 = m , the Bondi mass. The Komar 
mass of spinning Schwarzschild M^^ r = m + r^/r 3 has both global mass and a quasilocal 
term which falls off asymptotically. When m is set to zero, then the spinning 

Minkowski metric. 

The g sM metric has unusual properties. Unusual because there is quasilocal mass but 
no global mass, and there is an ergosphere but no event horizon. Values from the Komar 
superpotential provide global angular momentum, J sM = The quasilocal mass does 

fall off asymptotically. It arises from the kinetic energy of the rotating string fluid. As seen 
from conformal infinity, g sM is a large spinning object with no total weight. 

In this first work we focused on asymptotically flat seed metrics, and spun them up 
in a prescribed manner. Future work will allow other configurations, such as cylindrical 
symmetry, or extended matter, and the Bondi-Sachs metric form can be relaxed. 
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APPENDIX A: NULL TETRAD 
1. Spinning Schwarzschild 
The null tetrad which spans metric g^ CH of Eq. pjl is, with Q = r^/r 3 and A 2 = l — 2rrio/r, 
l^dx? = du, P«9 M = d r , (Al) 



n^dx" = {A 2 /2)du + dr, = d u - (A 2 /2)d r + Qd ( 

% sin t I % 

m^dx^ = ( -= — )du — {—=){d , d + ismfi dip), m^d^ = —^(d$ + - — ^d v ). 

\J2 \J2 y2r smv 

The Weyl tensor components are 

= 



3i Q sin $ 



2y/2 r 

q 2 = -I^-2(fiW# + i~cos#), (A2) 
r A 2 r 



2 r 



^4 = 0. 

The Ricci components are 



$00 = $01 = $12 = $22 = 0, 

$ 02 = ^ 2 sin 2 tf, 
o 

$11 = --$02, 

6A = TZ/A = $02- 



2. Spinning Minkowksi 

The massless spinning metric, with Q = r^/r 3 
g s ^dx a dx? = (1 - tt 2 r 2 sm 2 $)du 2 + 2dudr + 2VLr 2 sm 2 d dud^ - r 2 (d$ 2 + sin 2 $dy? 2 ) 
is expanded in a null tetrad as 

flS = l a np + njp - m a fhp - rh a m p (A3) 



13 



where 

l a dx a = du, l a d a = d r , (A4) 

n a dx a = du/2 + dr, n a d a = d u - (l/2)d r + Qd v , 

IT I Z 

m a dx a = —=(Qr sin $)du =(d$ + isnr&dip), m a d a = -= (d* H -du,) 

V2 V2 V2r sm-d v> 

This tetrad has 6 zero-valued spin coefficients: K = e = a = X = u = '-f = 0. The others are 

1 1 

P = — , V= -7T> 

r 2r 

n = ( — ^Ifisini?, t = —Ti, 

/"is geodesic since k = 0, and n a is geodesic since u = 0. 
The Weyl tensor components are 

*o = 0, 

. 3i s Q sin $ 



2sf2 r 



tf 2 = -hn 2 sm 2 ^ + i-cos^), (A5) 
^ 3 = -*i/2, 
^ 4 = 0. 

If there was a Bondi mass, it would appear in \& 2 at 0(l/r 3 ) since r is a valid luminosity 
distance. The Einstein tensor null tetrad components are 

G% = -mi{l a np + njp) - 2$ 02 (m a m /3 + m a m p ) + (2$ n - K/4)g$ (A6) 

where 

$00 = $01 = $12 = $22 = 



9„o . o „ C/ 

02 



$ n2 = -ft 2 sin 2 $ 



1 
2 



2 

$ n = --$02, 



6A = K/4 = $02- 

One can transform from the null tetrad to the locally non- rotating tetrad in Eq. ffTBj) : 

v — l/2+n, r — — 1/2 + n, ■& — — (m +m)/V2, (p = i(m — m)/v2. (A7) 
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APPENDIX B: INVARIANTS 



In general, the quadratic Riemann invariant (Kretschmann scalar) is related to the Weyl 
and Ricci invariants by 

Raf3fj,uR al3fJ/l/ = C a /3^C af3fJ/U + 2R a /3R a(3 — 7£ 2 /3. 
The quadratic invariants for g sSCH are, with Q = r^/r 3 

R% C X£n = + 24^(3fisintf) 2 + i^ST , (Bl) 

C$£C°g& = 48^ + 24^(3fisin^) 2 + ^Cr, 
R% CH R:£ck = ^ sin 4 K 2 sch = siu 4 , 
and the quadratic invariants for g sM are 

R B X„RT = -3(^) 2 + 2(^^) 2 + ^Vsin 4 0, (B2) 
C*^C$r = -3(^) 2 + 2(^^) 2 + 108^sin^, 

dsM r>a/3 243 4 . 4 2 ^^ri4 • 4 q 

R a pR S M = —V sm 0, K sM = —Q sm 0. 
For comparison, the Kerr invariant is 



pkerr -notfi^v ^kcrr f~iO.fi {iv 

■^afi/iu-^kerr ^ afi^v ^kerr 



Ar, m o l-a 2 cos 2 ,„ _ , 2n 4 4 

48 7^ 5 5-5nb(! - 14a 2 cos 2 + a 4 cos 4 

r 6 (1 + a 2 cos 2 0) 6V ' 



where a = a/r. 



APPENDIX C: POINCARE INEQUALITY 

An inequality due to Poincare jQ] shows that ft is reasonable. Consider the equilibrium 
of a rotating, self-gravitating fluid obeying Newtonian gravity. 

V 2 = AirGp 

for uniform density p. The fluid is spinning with constant Q about the z axis. The force on 
the liquid is 

F = -V[4 > -(n 2 /2)(x 2 + y 2 )]. 
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For equilibrium, the total force on the boundary of fluid volume V with outward normal n 
must be 



F-nd 2 x < 



V[0- (n 2 /2)(x 2 + y 2 )] -nd 2 x < 
J V 2 [<j)-(n 2 /2)(x 2 + y 2 )}d 3 x < 
J (4vrGp - 2Q 2 ) d 3 x > 0. 



The mass, volume, and average density are 

M = J pd 3 x, V = J d 3 x, p = M/V 

and so 

Q 2 < 2nGp. (CI) 

It is simple to verify that Q(r) = r^/r 3 satisfies Poincare's inequality for r > r$ and quasilocal 
mass M = r^/r 3 in Eq. (|23j) . The average density is p = M/(47rr 3 /3) = 3^/47? and Q 2 = 
r^/r 6 . 

APPENDIX D: BONDI-SACHS METRIC 



Asymptotically flat systems are described by the Bondi-Sachs metric 



9™dx*dx v 



Ve 



2h 



du 2 + 2e 2b dudr - r 2 H AB (dx A - U A du)(dx B - U B du), (Dl) 



B ttB, 



where outgoing null hypersurfaces are labeled bv u and coordinates are (x° = u, x l 



r, x 2 = x 3 = if). The Bondi-Sachs metric jl5( extends Bondi's original metric Q| to 
include ip dependence and has six independent functions {V, b, IT 9 , U^, y, q}, of (u, r, if). 
The rays of each u = const null surface are null geodesies x a (r) with tangent dx a /dr where 
x 1 = r is a luminosity parameter. Coordinates (p) are constant along each ray. The 
luminosity parameter is defined by r 4 sm 2 {} = det(^As) = det(r 2 HAB) where 



H 



AB 



e 2y cosh(2g) sinh(2g)sin$ 
sinh(2g) sin $ e~ 2y cosh(2g) sin 2 $ 
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The boundary conditions on the metric functions in the limit of future null infinity 1 + are 

rU A ^0, 6^0, y^O, q^O, V/r -> 1. 

Notation was chosen to avoid confusion between Sachs metric functions and Newman- 
Penrose spin coefficients: 

2y = 7 + 5 (Sachs), 2q = 7 — 5 (Sachs), b — (5 (Sachs). 

The u = const hypersurfaces have null geodesic tangent l a d a which is also hypersurface 
orthogonal as l a dx a = du. The twist of l a is zero and its expansion and shear are given by 

p = -e~ 2b /r, a = -e- 2b [{d r y) cosh(2g) + i{d r q)\ (D2) 

which follows from r as a luminosity parameter, and where the phase of a is determined by 
the choice of tetrad orientation. The Bondi mass aspect M(u, ip) is 

-2M = f» + $» + a u (A ), (D3) 

where the zero superscript indicates the leading coefficient in a 1/r expansion. The Bondi 
mass is the two-surface integral of the mass aspect over a topological two-sphere at X + 

M Bo ndi = j M + V°2 + d u (a°a )} y/=jjd&d<p. (D4) 

s 2 



[1] E.T. Newman and A.I. Janis, J. Math. Phys. 6, 915 (1965). 

[2] E.T. Newman, T.E. Couch, K. Chinnapared, A. Exton, A. Prakash, and R. Torrence, J. Math. 

Phys. 6, 918 (1965). 
[3] E.T. Newman, J. Math. Phys. 14, 774 (1973). 
[4] S.P Drake and R. Turalla, Class. Quantum Grav. 14, 1883 (1997). 

[5] J. Ehlers, in Les Theories Relativistes de la Gravitation, Eds. A. Lichnerowicz, M.A. Tonnelat 

(Colloques Internationaux C.N.R.S. No. 91, Paris 1962). 
[6] R. Geroch, J. Math. Phys. 12, 918 (1971). 
[7] R. Geroch, J. Math. Phys. 13, 394 (1972). 
[8] J. Bicak and A. Pravdova, J. Math. Phys. 34, 6011 (1998). 



17 



[9] E.N. Glass and M.G. Naber, J. Math. Phys. 35, 4178 (1994). 

[10] J.D. Lawrence, A Catalog of Special Plane Curves, (Dover, New York, 1999). 

[11] M.P. do Carmo, Differential Geometry of Curves and Surfaces, (Prentice-Hall, New Jersey, 
1976). 

[12] J.L. Synge, Relativity: The General Theory, (North-Holland, Amsterdam, 1964). 

[13] P. Szekeres, J. Math. Phys. 6, 1387 (1965). 

[14] H. Poincare, Acta Math. 7, 259 (1885). 

[15] R.K. Sachs, Proc. Roy. Soc. (London) A270, 103 (1962). 

[16] H. Bondi, M. van der Burg, and A. Metzner, Proc. Roy. Soc. (London) A269, 21 (1962). 



